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1.  INTRODUCTION 

A  single -determinant  form  for  wave  functions  has  been  used  extensively 
in  the  quantum  mechanics  of  many -electron  systems. 

In  the  usual  Hartree-Fock  method,  single -determinant  wave  functions 

are  used  with  a  restriction  that  two  electrons  with  spins  a  and  (3  are  put 

into  the  same  space  orbital.  A  single -determinant  wave  function  with  this 

2 

restriction  is  an  eigenfunction  of  the  total  spin  operator  S  .  The  Hartree- 
-Fock  method  has  proved  to  be  very  useful  in  the  theory  of  atoms,  molecules 
and  solids. 

In  the  unrestricted  Hartree-Fock  method,  single -determinant  wave 
functions  without  the  above  mentioned  restriction  are  used.  This  allows  us  to 
treat,  in  a  compact  form,  the  exchange  polarization  and  the  correlation  of 
electrons  with  antiparallel  spins.  However,  the  unrestricted  Hartree-Fock 
method  has  a  disadvantage.  Wave  functions  used  in  the  method  are  in  general 
not  eigenfunctions  of  the  total  spin  but  are  linear  combinations  of  eigenfunctions 
which  have  different  eigenvalues. 

It  is  of  some  interest  to  see  how  much  of  each  spin  eigenfunction  is 
contained  in  the  single -determinant  wave  function. 


2.  FORMULATION  OF  THE  PROBLEM 

Any  antisymmetric  wave  function  of  an  N-electron  system  can  be  written 
in  the  following  form, 

=  ,Im)0(Oi,O2,-.ON))  (2-1) 

where  ^  stands  for  the  space  and  spin  coordinates  fj  and  <r.  respectively 
of  the  i-th  electron.  A-  is  the  idempotent  anti  symmetrizing  operator.  By  using 
permutation  operators  P  and  their  parities  E  p  ,  the  anti  symmetrizing 
operator  A  is  expressed  as 

^  =|\jl  p  £p  ^ 

The  wave  function  $  is  in  general  a  linear  combination  of  pure  spin 

states: 


S.M 


$  =  15 


(2.3) 
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where  ^  is  an  eigenfunction  of 
and  M  respectively  (h  =  l). 


S  and  Si 


with  the  eigenvalues  S(S  +  1) 


This  decomposition  is  of  physical  importance  if  the  Hamiltonian  does  not 
involve  spin  operators.  An  expectation  value  of  a  spin-free  operator  f  can  be 
expressed  as 

<35,  S> 

_  r  (2.4) 

=sIMU)s'M  <f>S'M  ' 


where 


and 


CO 


S,M 


JlsvMLiO,  (zusm-D 


S,M 


(2.5) 


<$SM,f  ^s.M^ 


(2.6) 


when  <  ^S,  M  ’  ^S.M 
-free,  we  have 


>  is  not  zero.  When  the  Hamiltonian  of  the  system  is  spin- 


E  —  £  Es,m  '  (2«7) 

S,M 

This  equation  shows  that  at  least  one  of  the  energy  expectation  values  Eg  ^  is 
lower  than  E  unless  all  Eg  ^  are  equal  to  E  .  By  selecting  from  <I>  a 
suitable  spin  component  $g  ^  ,  we  have  a  wave  function  which  is  not  only  a 
spin  eigenfunction  but  which  has  a  lower  energy  expectation  value. 

The  analysis  is  also  useful  in  interpreting  the  function  $  ,  (2.  l).  This 

is  in  some  cases  (e.g.  a  single  Slater  determinant)  much  easier  to  handle  than 
its  components  ^  .  If  we  know  the  values  of  cjg  ^  and  ^  ^  M  '  we 

can  better  interpret  the  nature  of  the  simple  form  (2.  1).  One  of  the  basic 
problems  here  is  to  determine  the  weight  o)g  ^  . 

In  the  following  we  shall  investigate  a  special  case  in  which  the  spin 
part  ©  of  the  wave  function  (2.  1)  is  a  simple  product  of  spin  functions  y\(o\) 
which  are  either  a(cn)  or  p(cr.)  : 


.e(o1(a2.-,oN  )=y1(o1  )Y2(°2  )  ••  Yn(°n  )  • 


(2.8) 
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The  spin  function  0  is  an  eigenfunction  of  S2  •  If  we  denote  the  number  of  a 
and  (3  functions  in  (2*8)  by  Nq  and  respectively,  the  eigenvalue  of  S  Z 
(M)  is  expressed  as 

M  =  -j(Na-  Np)  .  (2-9) 

In  order  to  evaluate  o>g  ^  ,  it  is  convenient  to  introduce  the  spin 

operator  Og  ^  which  projects  out  the  component  of  the  pure  spin  state: 


fs2-s(s+i)}cv  0 . 

{SK-  M  }  Os,m  =  0  .  (2. 10) 


We  note  that  this  operator  works  only  on  the  spin  part  of  a  wave  function.  For 
any  wave  function 


$  =  ICj  Vj(i,,r2.  "•  ,rN)0i(or,o2.-.oN  )  , 

i 

w  <£  is  expressed  as 
0,  M 

os,Ma>=  icivi(os.Mei) . 

Since  the  wave  function  under  consideration  is  an  eigenfunction  of 
with  the  eigenvalue  M  ,  we  shall  drop  the  subscript  M  in  the  following  unless 
it  causes  some  ambiguity. 

The  weight  u>g  for  the  wave  function  $  can  be  written  as 

<os  $,o5  $>_<mos4'i/0> 

Ws  "  <  $  y  $)  =  (^0,4^0) 

i  Ep^,  px  vxe,0s  P°e> 

P  _ 

1  e^p^xe/e)  .  (2,11) 

P 

x  cr 

where  P  and  P  denote  the  corresponding  permutations  of  the  space  and  the 
spin  coordinates  respectively.  Similarly  we  obtain  the  expectation  value  of  a 
spin-free  operator  f  , 
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<os  $,f  os$>  <^afos4^e> 

/fv  -  — £_J - £ — - £ - 

s  .<os  $,os  $>  (mosAii/e> 

i  £P  w  p^xaOsP0©) 

i  EP(^px^>  (e,os  p°e>  •  (2,12) 

We  shall  derive  an  explicit  formula  for  <0,  OgP^B  >  in  the  next  section. 
In  sections  4  and  5,  the  weights  u>g  derived  from  a  single  Slater  determinant 
are  given,  and  the  behaviour  of  the  Wg -values  for  large  N  is  discussed. 

3.  CALCULATIONS  OF  <  0,  OgP0”©  > 

Because  we  assumed  the  form  (2.8)  for  the  spin  part  ©  ,  for  any  permu¬ 
tation  P  we  can  find  a  permutation  Q(P)  which  brings  ©  and  P0"©  into  the 
following  forms: 

Q°0=aO)a(2)  -  a(n+M)p(n+M+1)  •p(2n)  s  f0  ,  (3- 1) 

Q°P°©=  a(1)- a(n+M-i)(3(n+M-i+1)~-  p(n+M)x 
xa(n+M+l)---a(n+M+i)p(n+M+i+l)"'(3(2n) 

(3.2) 

5  f  i  .  (  2  n  =  N  )  . 


Here  the  integer  i(P)  is  the  number  of  a -functions  in  ©  which  are  changed 
to  (3 -functions  in  Pa©  .  The  number  i(P)  is  uniquely  determined  by  the 
given  permutation  P  • 

Using  the  commutability  of  Og  and  Qa  ,  we  obtain 

<0,OsP°©>=  <Qo0,Q°OsP°©> 

=  (Q°  0,OSQ°P°0) 

=  < fo,Osfj)  =  cs  j 


(3.3) 
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In  order  to  calculate  the  value  cc  .  ,  it  is  convenient  to  divide  the 

o  ^  x 

total  number  of  electrons  into  four  groups  A,  B,  C  and  D  .  A,  B,  C  and  D 
stand  for  the  first  n  +  M  -  i  electrons,  the  second  i  electrons,  the  third 
i  electrons  and  the  last  n  +  M  -  i  electrons,  respectively; 

A  1,  2,  . • . ,  n  +  M  -  i  , 

B  *n+M-i  +  1,  n  +  M-  i+2,  ...»  n+M, 

C  n+M  +  1,  n+M+2,  n+M+i, 

D  n+M+i+1,  n+M+i+2,  2n. 

For  example  when  we  write  CD,  this  means  the  combined  groups  of  C  and 
D  ,  i.e.  the  last  n  -  M  electrons. 

We  introduce  the  symbol  Y^(s,  m;  p )  for  one  element  of  an  orthonormal 
complete  set  of  simultaneous  spin  eigenfunctions  of  S  and  5^  with  the  eigen¬ 
values  s ( s  +  1)  and  m  respectively#  K  denotes  an  electron  system  which  can 
be  any  of  the  groups  defined  above.  \x  specifies  one  of  the  spin  functions  with 
common  s  and  m  ,  to  differentiate  degenerate  functions.  We  adopt  the  usual 
convention  for  the  relative  phase  of  these  functions 

( Sx±  i  SY)  YK(s/rn;|i)=7(s?m)(s±m+1)YK(s,m±1;jJL) .  (3.4) 

c2 

As  the  element  of  the  set  which  has  the  highest  eigenvalue  of  3 
(s  =  Nj^/2,  Nk  is  the  number  of  electrons  in  the  group  K  ),  we  choose  the 
following  function: 


VK(  Nk/2,Nk/2)  =  aa  •••  a  .  (3.5) 

We  may  drop  \i  in  this  case  since  (3.  5)  is  the  only  function  in  the  set  with  the 
eigenvalues  s  =  m  =  N^/2  .  From  (3.4)  and  (3.5),  we  obtain 

Yk(Nk/2,-Nk/2)=3P  -  p  .  (3.6) 

An  orthonormal  set  of  spin  functions  of  the  system  ABCD  can  be  obtained 
by  coupling  YAB  and  Y^  in  the  following  manner: 

Yabcd( s' m:  M- )  =  VABcD(s«m:  s' SV  ^ 

^b^s1/  rn'*»  fi  )'^D(s/rri;  3! )  x 

x  (s'm '  s'!m"  I  s',  s ,  s,  m  ) .  ^3> 
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Here  (s!m',s','m"  |  s',s','s,m)  is  the  vector  coupling  coefficient 

(s;m‘,  s"m"l  s',s*,'s,m)  =  6(m’+m"-m)* 

f  (2s+1)(s'+  s'- s)! (s’- m')! (s'- m")!(s+  m)!(s-  m)(j1//2 
„(s'+  s+  s+1)!(s-s+s)!(-s+s+s)!(s,+  m’)!(s"+  nrf)! 


*  I  (-0 

V 


V  +  s-  m'  (s'+  mV  v  )!(s+  s  -  m-  v  )! 


v!(s-  m'-v)!(s- m-v  )!(s- s+  m'+  v)!  1  (3.8) 


Since 


=  £  ^abcd(S,  ))  (Vabcd(5,M;(j.  )  t  (3.9) 


we  obtain 


cs, i  -  £  ^0 -yabcd(S,  M;(i  )>(  YABCD(S,M;  pi  )(fj>. 


(3.10) 


The  functions  £q  and  are  expressed  as 


f  _  v  (  n  +  M  n+M  )  v  (n~M  -  n-M  n  •  ,3  n) 

To  -  yab^  2  ’  2  '  Ycd^  2  •  2  '  -  ( 


and 


x  _  v  (  n+M-i  n+M-h  v  /  i  i  \ 
7i  ~  V  2  .  2  B  k~27~2 

«*<*.*)  vfatn^Mri -nzM=i ) , 


(3.12) 


Therefore,  the  terms  of  the  right  side  of  (3.  10)  vanish  except  for  the  following 
Y  . 


ABCD  * 


V 

YABCD 


YflR(— ■  m')Yrn(^  m"  )■ 


2  .  2 


m'nV 


2  .  , 


'  2 


■  2 


(3.13) 
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(3.  14) 


W11^, m )  n£A(r>*?  nV)  V }> 


*(  ~y~i. m" I y.-P — y — ^.yy-.m)/3  5> 


Using  (3.8),  (3.11),  (3.12),  (3.13),  (3.  14)  and  (3.  15),  we  find 

9sj*  <fo.YABcD(S.M:^'^)><YABCO<S.M:^,-^).fi> 


,n  +  M  n  +  M  n-M  nrM  ,  n  +  M_  rv-_M  ^  x„ 
25  \  2  '  2  '  2  2  2  *  2  /  ^ /  ’v»  / 


2  '  2  -  2  • 
n+M-i  n+M-i 


2  -  2 


n+M-i  _i_  n+M  n+  M-  2 
2  <2<2.  2 


, i_  n-M-i  n-M-i  . j_  n-M-i  n-M  n -M-2 i  A)C 
*4  -  2  <  2  <~  2  2 «  2-2-  2  j 

/  n  +  M  n  +  M  -  2i  n-M  n  -  M  -i2  i  .  n+M  n-M  <-  m ) 

*  \  q,  o  ,  0/  0  9'  ?  *  ^  i 1  ' 


_/0c  ,x(n-M-i)!i!  J-v  (n  +  M  -  i  +  v  )!(S- M  +  i-v  )i 
+  (n  +  S  +  1)!v  v!(i-v)!(S-M-v)!(n- S -i  + v  )! 


■■(oci^n-M-D.^iML  ny  {(S-M+v  )!}  _  (3  17) 

-(2b+l)  (S_ M)|  ?  V (S-M+v  -i)!(n-S-v)!(2S+1  +  v)!  .  (  1 

It  is  convenient  to  use  the  expression  (3.  16)  in  calculating  the  values  Cg  ^ 
for  some  special  cases.  For  example, 
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4.  APPLICATION  TO  A  SINGLE -DETERMINANT  WAVE  FUNCTION 

When  a  wave  function  $  is  a  single -determinant  wave  function,  the  space 
part  ¥(I|  IT2  )  is  a  product  of  one -electron  functions: 

,iN)  =  ^  (r t  )ip2(r2)  •••  ^N(iN)  .  (4-^ 

We  denote  orbitals  associated  with  a -spins  by  <^,  <j>£,  •  •  •  ,  ^+^4  anc^ 
those  with  (3 -spins  by  cp^,  cp^»  •  •  •»  9n  ^  •  Without  changing  the  total  wave 

function  $>  ,  we  can  transform  the  orbitals  {<)>}  and  {qp}  to  (4>'}  and  {cp 1 } 

3) 

so  that  the  only  overlap  remaining  is  between  the  pairs  <J>/  and  cp/  * 

n+M 

<j> /  =  Z 8 ..  ,  i  =  1,2(-*,n  +  M  # 

Y‘  j=1  u  J 

n-M 

<p.'  =1  b..ip.  ,  i  =  , 

.  1  j-i  'i  J 

;  (<!>;, 4>/>  =<*;.*/>  =  6^, 

«m  a;.  m°- 

Using  these  transformed  orbitals  for  the  space  part  ¥  ,  we  find  that 
the  inner  product  <  Y,  Px¥  >  vanishes  except  when  PX  is  a  product  of  some 
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interchanges  of  the  pairs  <p/}  •  When  P  interchanges  the  electrons  of 

t  pairs  <p'kl)  -  {<l>'k2>  ^kZ  ^ .  {<t>'kt  >  <Pkt}  •  we  see  that 

(qj,PK\V)  =  Xkl  Xk2  •••  Xkt  , 

Ep  =  (-1  )!  (*•  3) 

i(P)  =  t  . 

It  follows  from  (2.11),  (3.3)  and  (4.  3)  that 

n-M  ,  ,k 

“s  =  I  (-DAk  cs  k  , 

k.O 

where  is  defined  by  the  coefficients  of  the  polynomial 

n-M  n-M  l 

n  (i  - x k x )  =  I  (-1)Akxk  . 

k  »  1  k  =  0 

For  example, 

A-o  =  1  - 

Ai  =  I  X  k  , 

a2  =i(r  \f-i  a,  =a1a2+xix3+  •■■+x,xn.M+ 

+  x2x3  +  ■  +  Xn.M.1  Xn_M  , 

An-M~  ^  • 

k 

In  order  to  calculate  (4.4),  it  is  convenient  to  introduce  the  following 
polynomial  W(x) 

,  %  n-M  n-M  l  l 

W(x)  =  TT  (1  -  (1  -Xk)x  }=  I  (-1)  Bk  x  . 

k»1  k=  0 


(4.4) 


(4.5) 


(4.6) 
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By  comparing  (4.  5)  and  (4.6),  it  is  found  that 


a  =  v  (_  i  y _ Ln _ m _ ih _  g . 

k  >0  (n-M  -  k)!(  k-  j)!  bj  ’ 


(4.7) 


The  weights  ojg  of  a  single -determinant  wave  function  are,  therefore,  expressed 
as 

.  l  >  1  J+k  (n-M  - j) !  0 

k  =  0  j  =  0  (n-M  -  k)!  (  k  -  j )!  J  s,k 

/  o c  .  1)  nyS nrM  y‘M+V/  1  )V  +J  +k 

-(2S+1)  (S-M  )!  v=o  ifo  k  =  i  ~ 


{(  S  -  M  +  v  )!j(n  -  M  -  j) 


Bi 


(n-S-v)!(2S+1+v)!v!(k-j )!(S -M  +  v  -k)!  1 


.vV  {(S-M  4V)!)2 

~Ub+i;(stT^!v,o  j(2S  +  1  +  v)!  v!  0S-M*v 


(4.8) 


This  can  be  written  as  an  integral, 


,,  /  i\S-m  2S+1  f  Ys-MM  ^s+Mds'MW(x)  , 

(S4«i  (  ’  '■ dX  ’  (4'9) 


When  W(x)  is  given,  we  can  derive  using  the  above  equation  (4.9). 


*) 


Since  the  first  order  density  matrix  completely  determines  the  original 

single -determinant  wave  function,  it  should  also  determine  the  weight 

u>  .  We  derive  an  explicit  expression  for  W(x)  in  terms  of  the  first 
b 

order  density  matrix,  Wg  can  be  derived  from  W(x)  by  using  eq.  (4.  9)* 

The  first  order  density  matrix  of  the  wave  function  under  consideration 
has  the  following  form: 

p(I,  if)  =  p+(r,Jr')a(o)  a(o')  +  p_(i.,r')  p(o)  (B(o'), 
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where 


Since 


n+M 

p+  =  Z'  4>i>  <4>i 


n  -M 

p.  =.i  «Pi><q>;  . . 


<ipi\  {(1  - x ) p.  +  xp.p+}  qDj'>  =  6jj{1  -(1  -Xj)x) , 


W(x)  can  be  expressed  as  an  expectation  value  of  an  n-M  particle 
operator  K(x): 


W  (  x)  =  (W,  K(x)  40  . 

where 

K(x:i1>f2.-.In.M:Il.l2.-.In4<)  = 

n-M 

=  TT  {( 1 - x)p.( Ij , f  j )  +  X  p.p+  (ij  ,  ffj*  ) } 
i  =  1 


and 

V  "(*1  .  I2;  •  •  )  =  /( n  )'^i  ( *1  XpK r2)  •  ip’n-  • 

It  is  seen  that  A  K(0)  is  a  projection  operator, 

A  K  ( 0)  =  40  <  r  , 

and  K(0)  K(x)  =  K(x)  .  Therefore,  we  obtain  an  expression  of  W(x)  in 
terms  of  p+  and  p  only: 

W(x)=<r,K(x)U0  =tr4K(0)K(x)  =  tr4K(x)  . 
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5.  BEHAVIOUR  OF  FOR  LARGE  N 
We  discuss  the  case  when  all  X^'s  are  equal  to  X  .  Then 

W(x)  =  fl  -  ( 1  -  \)x)n  M. 


(5.1) 


Putting  (5.  1)  into  (4.  9)»  we  obtain 
(1) 


{D,  i ;  into  we  outam 

M- (2S+1  ( 1  -xri1-°  -X)xr 


’dx. 

(5.2) 


The  asymptotic  form  of  (5.2)  is  derived  in  the  appendix  when  S  -  M 
2/3 

is  small  compared  with  N  '  , 


^  ~  (1.2x^n-M'S~f*'e-P  ' "^<S  '  "  _1  XS'M) 

1  -  X 

_  2M+(1„rX)2(n-M)i.  (S-M  +  1  )(S-M  +  2)- 
2  ( 2M+(1  ~\)(n  -M)r 

-1m7(T^m)(s-m)(s-m+1))-  <-> 


The  equation  (5.  3)  may  be  regarded  valid  over  all  possible  values  of  S  ,  since 

both  left  and  right  hand  sides  of  (5.3)  decrease  rapidly  as  S  -  M  becomes 

l/2 

large  compared  with  N  '  . 


5  a)  Orthogonal  case 

Putting  X  =  0  in  the  expression  (5.2),  we  obtain 


(5.4) 


It  is  interesting  to  note  that  u>s(0)  is  proportional  to  the  number,  fn  g  »  °* 

linearly  independent  spin  functions  for  given  n  (  =  N/2)  and  S  , 


WsCH.-ir-ki 


I  fn.  S' 
S'=  M 


(5.  5) 
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The  asymptotic,  form  of  coq(0)  is,  from  (5.  3), 


“S(0)"(2S+1)(57^'exp{’2(^(s-M-1)(s-M)- 


1 


2(n+  M ) 


(S-M+1)(S-M+2)}  . 


(5.6) 


For  M  =  kn(0  <  k  <  l)  ,  •  cog  can  be  approximated  by  a  geometrical 
sequence, 


coc  ~  (2S+1) 


(n-M) 


S'M  2k/1-k's"M 


(n  +M )5_M+1  1  +  k  \1  +  k 


(5.7) 


It  should  be  noted  that  (5.7)  does  not  contain  the  number  of  electrons  explicitly. 
Therefore,  for  a  fixed  k  ,  o>g  does  not  change  much  when  N  increases. 

For  M  =  0  , 


4  (2S+1 )  /S+S+1 

ws  -  — — exp  - 


n 


n 


(5.8) 


2  3 

From  (5.8),  the  expectation  values  of  S  ,  S  and  S  can  be  calculated  as 


2  2 

<S>  =  I  5  ws~  J-^exp(-~ -)dS  =  Jrm/2  . 
<S2>  =  Z  S2^-  J^exp(--^2)dS  =  n  , 

<S*>  =  I  S3u^~  exp(-  -^-)dS  =  3 /n3n / 4 


t) 


The  exact  values  are  as  follows, 

<s>  - _l 

2  (2  n) ! 

<S2>  =  n  +  1  -<S>  . 


<c3>  _  1  .  22n~2(6n  -5)  n  (ni)2 
4  (2n  +  1)! 
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In  figures  1  and  2,  Wg  is  plotted  as  a  function  of  S  and  S/n 

respectively  for  N  =  10,  100  and  1000,  using  (5.8).  We  see  from  these 

figures  that  when  N  increases,  spreads  towards  bigger  N  .  However, 

for  large  N  ,  appreciable  weight  o>  is  localized  around  the  value  S  ~ 

o  max 

~  0.  5*VN  •  The  second  moment  of  the  distribution  around  the  average  <  S  > 
is 

<S2>-  <S>2  -  (1  ~)n  ~  0.107  N  . 

4 

5  b)  Non -orthogonal  case 

For  M  /  0  , 

u)s(n,M,X)  ~  u>s((1  -X)(n-M)  +  M,  M,0  }  ,  (s-9) 

since  the  factor  (2S  +  l)(n  -  M)^  ^/(l  -  X)(n  -  M  +  2M/ 1  -  X)^  decreases 

rapidly  compared  with  the  exponential  part  in  (5.3),  which  may  therefore  be 
regarded  as  1  . 

For  M  =  0  , 

^nM|-Mne4S11-l)n' 

~  (i)g  { ( 1  -X) n, 0 , 0  }  .  (5-io) 


It  follows  from  (5.  9)  and  (5.  10)  that  the  weight  in  the  non -orthogonal  case  for  N 
electrons  can  be  approximated  by  the  weight  for  2  { ( 1  -  \)(n  -  M)  +  M}  electrons 
in  the  orthogonal  case: 

WgCn.M.X)  ~  o)s{(1-X)(n -M)  +  M,M,0}  (5.11) 


When  the  overlap  X  increases,  the  distribution  shrinks.  This  can  be 
seen  from  figures  3  and  4,  where  for  some  X  values  are  plotted  as  a 
function  of  S  for  M  =  0  and  M  =  0,2  n  respectively. 
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APPENDIX 


We  define  a  function  I  by  the  integral, 

Ia,p,y  (z)  =  — j  fxa(1-xf(1-ZX)Ydx 

Ct!  Jq 


(A.  1) 


Then  the  weight  (5.  2)  can  be  expressed  as 


vS-M 


■M*)-  ttSd )(-Q7(^:Lsy—  ■s-M.s.M.n-st1  -x)-  (A-2) 

The  maximum  of  the  integrand  in  (A.  1)  is  given  by  one  of  the  roots  of 
the  following  equation, 


a.  _  P  zy  =  0 

x0  1-x0  1-zxq 


(A.  3) 


Therefore, 


„  (1  +z)a  +  p+ZY-i((1+z)a+  (3  +zY)2-4az(a+p+y)  ,, 

X° - 2Z(ai?+7)  -(A-4) 

The  order  of  xq  is  the  same  as  a/(3  +  zy). 

Expanding  (l  -  x)^(  1  -  zxj^exp  —  in  terms  of  x  -  x  ,  we  obtain 


1 


a! 


x  (1  -xf  (1  -zx)  = 


=  {^q+  a^Cx  -;xb)  +  a3^X  -X0^+  ^  ^*P  (  ocx/Xq). 


(A.  5) 


Here 


and 


a0  =d  -x0)p(1  -zx0)Yexp  a  . 

-  =  _  ±of  _j _ +  _ih _ ) 

'2  2  \(1  -  x0)2  (1-zx0)2/  , 

ak  =.  a0  0  (affl  x0"(^l)  . 


(A.  6) 


(A.  7) 
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Int  eg  rating  (A.  5)  over  the  range  (x  =  0,  x  =  oo)  ,  we  arrive  at  an 
asymptotic  expansion 

I  ~  3q  bo  4-  cl 2  (A.  8) 

where 

bk  =~  fxa(x-  x0)kexp(-ax/x0)dx  ,  (a. 9) 

a»  j0 

b0  =  (x0/a)a+1, 
b2  =  b0  Xo(a+2)/a2 , 

bk  =  b0O  (  Xq  ofnr)  )  .  (a.  ii) 


From  (A.  5),  (A.  7)  and  (A.  11),  we  find  the  order  of  magnitude  of  a^b^.  : 


ak  kk/a0  bo  = 


K 

0  (xJ  )  k  :  even  . 

0(  X0Ri//p+  Z  Y)  k  :  odd  . 


(A.  12) 


Remembering  that  the  order  of  xq  is  the  same  as  that  of  a/(p  +  zy)  and 
using  (A.  12),  we  can  take  the  first  two  terms  in  (A.  8)  in  order  to  calculate 
log  Wg  with  the  accuracy  of  order  (p+  Zy)  .  Then  we  obtain 

j  C  _  *o  (a+2)f  (3  .  z2y  l']  „ 


2a2  1(1  - x0) 2  (1  - zx0); 


y  U.T-  l 

x  1  -*o  f(  1  -zx0)Yexpa  . 


(A-  13) 


The  expression  (A,  13)  is  substituted  into  (A.  2),  and  after  some  manipulation, 
we  obtain  (5.  3). 
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